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This appendix to the paper by Werpachowska and Dietl Theory of spin waves in ferromagnetic 
(Ga,Mn)As^- is a response to inquiries about the derivation of the Lowdin calculus and its numerical 
implementation. It presents the detailed version of the Lowdin calculus for the multiband Hamil- 
tonian taking into account both the p-d and s-d exchange couplings. This should explain how to 
implement the sp-d exchange interaction in the tight-binding computational scheme. I also include 
the complete procedure of the Bogoliubov transform for systems with spatial inversion asymmetry. 

I. EFFECTIVE HAMILTONIAN 

We consider the ferromagnetic phase of a system consisting of P carriers and N magnetic lattice ions, described by 
the Hamiltonian Hq and coupled by the sp-d exchange interaction %' , 

p N 

H = Ho + H l =H + J2Y, S *- V( r * - R j) > (!) 
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o 

^ ■ where Sj and Sj are the i-th carrier's and j-th ion's spin operators, while and Rj are their respective positions. 
The matrix elements of / between the s- and p-type basis functions (a and /3, respectively) determine the strength of 
the p-d and s-d exchange interaction, respectively. They arc defined as 
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a = [ d 3 rS(r)*J(r - R)S(r) , 
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[ d 3 rX(r)*/(r-R)I(r)= f d 3 rF(r)*/(r — R)F(r) = f d 3 rZ(r)*I(r - R)Z(r) , (2) 
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, where S(r), X(r),Y(r) and Z(r) are the periodic parts of the s- an p-type band wavefunctions at k = 0. Their values 

for (Ga,Mn)As are N /3 = —1.2 eV and N a = 0.2 eV, where N is the concentration of cation sites^— . In the absence 
CS| ■ of external fields, Ho depends on the carriers' degrees of freedom only. 

The different signs of p-d and s-d exchange integrals lead to antiferromagnetic coupling of the ion spins with the 
valence band p-type holes and conduction band s-type electrons, respectively. Within the mean-field and virtual 
crystal approximations, these couplings produce a Zeeman-like spin splitting of the energy bands, described by the 
additional term in the one-particle Hamiltonian 



h = h + xN S(pP p + aP s )s z , 

where ho describes the host band structure, x is the fractional concentration of Mn ions with total spin S = 5/2 and 
P p , P s are projection operators associated with the p- or s-type band wavefunctions at k = 0, respectively, which 
commute with the spin operator. As the s-d interaction strength a is much smaller than /3, the latter determines the 
strength of the ion-carrier coupling. It is convenient to write h as 



H ' h = ho + As z P cx , (3) 

where A is the exchange spin splitting of the heavy hole p-type bands in the T point, equal -0.15 eV for x = 5% 
and rescaled linearly for different Mn concentrations, and P cx = P p + jP s commutes with the spin operators^ With 
the above definition of P ex , h is directly applicable to the k ■ p band structure calculation methods we have used in 
Ref. 0], as their basis functions are exactly the ones that operators P p and P s project on. However, the tight-binding 
approach uses localised, atomic-like orbitals wavefunctions as a basis. The Mn ions are coupled to p-type orbitals of 
As and s-type orbitals of Ga. The spin splittings applied to these orbitals must take into account their weights in the 
s- and p-type band wavefunctions close to the centre of the Brillouin zone. The operator P ex is therefore redefined to 

P — A -1 (p P As + p P Ga ) 

where A 5 % = —0.15 eV is the heavy hold band spin splitting for 5% Mn concentration, and Pp S and P s Ga are the 
projection operators on the respective orbitals, and the effective orbital spin-splittings for 5% Mn concentration are 
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Parametrisation e p [eV] e s [eV] 
Jancu TBA -0.2764 0.0548 
Pi Carlo TBA -0.1940 0.0356 

The values of the above constants were chosen so that the heavy-hole splitting A matches the one resulting from the 
k ■ p method (see Fig. 1 in Ref. [!). The weights of the projection operators are almost the same in both methods, 

&p, Jancu ^ &p, Di Carlo ^ P 
t-s, Jancu c-s, Di Carlo Q 

The dynamics of magnetic ions coupled to the system of carriers requires a self-consistent description, which takes 
into account how the carriers react to the ions' magnetization changes. Therefore, we use the Lowdin perturbation 
method specifically adapted for multiparticle Hamiltonians, 6-9 to derive an effective Hamiltonian T-L cS for ions only. 
The perturbation in question is the sp-d interaction between carriers and ions. 

We choose the multiparticle basis states of H as M ® T. The ion part M has spins quantized along the z direc- 
tion, while the carrier part T is a Slater determinant of P one-particle eigenstates ipKm °f the one-particle carrier 
Hamiltonian ([3]). The subscripts k and m denote the wave vector and the band number, respectively. 

The Lowdin calculus consists in dividing the multiparticle basis states into two subsets, A and B, 



n = 



Haa Hab 

/H.BA Hbb , 



Set A contains all states M (£) To, where Tq is a Slater determinant of the P lowest eigenstates of h. Set B contains 
all the remaining states, in which at least one carrier is excited above the Fermi level. We construct the effective 
Hamiltonian for the states from set A only, adding their coupling with set B as a second order perturbation, 

Kn' = (^o)nn' + H' nn > + V] "" - , (4) 

7rr„ tii — n n "n" 

where H nn ' — {n\%\n') (similarly for (Ho) n n' and H' nn i) and n,n' <E A. The term (Ho)nn' is independent of the ion 
configurations, so we set it to zero for simplicity. Thus, the effective Hamiltonian T-L cS depends only on the ion degrees 
of freedom but, thanks to the Lowdin method, takes carrier excitations into account and can be used to calculate the 
spin- wave dispersion in a self-consistent manner. 

The variational part of the presented method consists in searching for the energy E somewhere in the region for 
which we want accurate results, which in our case are the lowest eigenenergies of H, in particular the groundstate. 
For a known average spin splitting A, we can set E to the total energy of the carrier multiparticle state Tq, Er — 
S(k m)er -®km, where -Ekm is the eigenenergy of i/'km, and the sum goes over all occupied eigenstates in Tq. The states 
n" are of the form M"®T" . T" ^ Tq. To simplify the sum over n", we approximate the diagonal matrix element H n "n"> 
which depends on both M" and T", by the total energy of the multiparticle carrier state T", Ep" = m )er" ^km- 
It describes the interaction of T" with the average configuration of the ions' spins corresponding to the spin splitting 
A. We can thus write the Hamiltonian ^ in the following form: 



nf n , =H' nn , +J2 E Hnn " Hn " n ; , n,n> e A . C,] 
M" r"^r r ° r " 



The factor U n n"Un"n> under the sum can be written as (M<g>r |'H|r"<g>M")(M"<g)r"|-H|rn«)M'), where n = M<Z>r 
and n' = M'tgiTo. Since the denominator in Eq. ([5]) is independent of M", summing over M" is equivalent to inserting 
an identity operator in the ion Hilbert space, which allows us to write the last term as 



E 



(m ® r \n\T")(T"\n\T ® a/') 



r" # r Sr » - E r" 



We can thus treat T-L eS as a Hamiltonian acting on ion states only, 

- (M(g>r \', 

Et — Er 



U% M , = (M|«-|M<) = (M « m 9 M') + Y. <M8r °'T ><r ^' r ° 8M ' > ■ («) 

r"#r 



Since the sp-d exchange term in H, which produces the non-diagonal matrix element (M <S) r |%|r" ® M"), is the 
interaction of a single carrier with an ion, the only T" states which have a non-zero contribution to the sum over T" 
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in Eq. ([6]) are those which are created from Tq by just one excitation, -0km — > VVm' from below to above the Fermi 
level — we denote such states by Fk m k'm'- Hence, we have Ef — -^t , , — Pkm — -E/k'm' an d Hamiltonian © can 
be written as 

w$ M , = (M®r \H'\r ®M') + J2J2 %^^^( M ® r ol^|rkmk' m ')(rk m kw|H|r ®M') , ( 7 ) 

kk mm 

where /k m is the Fermi-Dirac distribution (this generalises the calculation to finite temperatures). 

To integrate out the carrier degrees of freedom in (0, we need to calculate the matrix elements (t/j^ m \s a I(r — 
R-jOIVVm'); where a = +,—,z and [s + ,s~] = s z by convention. Assuming that the function I(r — Rj) vanishes 
quickly outside the unit cell and that e lkr and e lk ' r vary slowly over the same range, we obtain 

(VwK/(r - RjMwrw) = e i( - W -^(u km \s a I(r - R^K'm') , 

Because of the spin-orbit coupling we cannot separate (wkm|s CT -f (r — Rj)l u k'm') into a product of the spin and spatial 
matrix elements. To overcome this problem, we write Uk m as J2 S ^km^s using the spinor basis ip s = |f), where 
u Lm i s a purely spatial wavefunction, and 

(u kTO | S ff /(r-R J 0|w k w)=^^ s |s CT |^)(< m |/(r-R J -)|<w) • 

S3' 

The operator P ex projects the wavefunctions on the basis functions which are coupled to the magnetic ion via the 
exchange interaction described by the function /. Hence, we can write that 

which leads to 

(«km|s°"-f(r-Rj)|ttkw) = cy^(^ a |s°'|^ s ')(Mkml- P exl^kw) = c («km|s°'Pcx|uk'm') , (g) 

ss' 

as P e x commutes with the spin operator. The scaling factor c can be derived from the condition that for the system 
fully polarised along the z axis and in the mean-field and virtual-crystal approximation, we have 

JV 

S 2j(ltkm|s Z f( r _ Rj)|«k'm') = A(u km |s 2: Pcx|uk' m '> , 
i=l 

which can be simply generalised to arbitrary temperature by replacing A and S by their temperature-dependent 
versions. Inserting ([5]l. we obtain 

iV5c(wkm|s 2 Pcx|uk'm') = A(w k m | S^Px | W k 'm' ) • 

But A = NS/3/V, where /3 is the p-d exchange integral defined in Eq. ©, so c = A/(NS) = (3/V and 

(^K^r-R^lVkw) = |e l ( k '~ k )- R ^ km k' m ' • 

where s^ mk , m , = (ukm|s CT P e x|uk'm')- The above expression is proportional to /?, but the operator P cx inside the 
matrix element takes into account the difference between p-d and s-d exchange integrals. 
We can now write Hamiltonian ([5]) using ion spin operators, 

JV JV N 

"•" EE H JSJ + EEE ntfsjsf . 0) 

a j=l era' j=l j' = l 

The coefficients HJ and H??, are given by 

R aa' _ ^ V V /k"^ 1 ~ /k'm') j(k'-k)-(R,-R,,) o- V (, jy, 

JJ — y2 / y / i £j _ £i e s kmk'm' s k'ffl'km ! V ±u ^ 

kk' mm' 
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where due to the condition T" 7^ Tq in Eq. j5), for k = k' the summation goes over m 7^ m! , and 



fr-£Wf, ^ _MV V - h,m') , a z a z , 

j y / j / j J™'kmkm y / , / , f _ jg^ l 5 kmkm ,; 'km , km T 'km'km^kmkmV ■ (11J 



To obtain the above expressions, we substituted the carrier-only part of Hamiltonian Hq ([TJ by the sum of P one-par- 
ticle Hamiltonians ho from Eq. ([3]). We also used the formula (V'kml^olV'k'm') = ^kk'(^mm'-^km _ ^ s kmk'm')' obtained 
from Hamiltonian h 

Using the Lowdin perturbation-variational calculus, we have thus described the problem as a lattice spin system 
coupled by the effective exchange interaction, integrating out the carrier degrees of freedom. The physics of the 
carriers is embedded in the effective Hamiltonian H cS (J9j) , and is responsible for the long-range nonlocal character of 
the mutual interactions between magnetic ions. 

The first term of H contains the operator responsible for the mean-field generated by the carriers acting on 
the lattice ions (the first part in ITTj). plus the correction arising from the inter-band transitions. The latter can be 
associated with the Bloembergen-Rowland exchange mechanism j 10 i 11 A small but worth noting result of the next 
section is that in the small oscillations approximation it cancels with the contribution from the second term. 

The second term of T-L cS describes the effective long-range exchange interaction between the lattice ions. The nature 
of the interaction, mediated by the carrier between two ions, is reflected in the appearance of the squared j3 constant. 
The s-d exchange interaction characterised by the a constant is incorporated in the s^ mk , m , matrix elements via the 
P cx operator. The fraction with the resonance denominator resulting from the perturbational approach dampens the 
influence of the distant energy bands, while the biggest contribution to the sum comes from the states in the vicinity 
of the Fermi level Ef- The denominator looks dangerous, as it may cause the fraction to diverge in the presence of 
the energy bands' crossings, which would make our perturbation calculus invalid. However, the effective Hamiltonian 
for ions depends on the average of these factors, and will be shown (at least in the small oscillations approximation) 
immune to this problem in Sec. |Hj 

Contrary to the original RKKY range function 

p(E F )k F J§ sm(2k F r) — 2fc F r cos(2fc F r) 
W = 2^ (2k F r) 4 

(where p{Ep) and k F are the density of states at the Fermi level and the Fermi wavevector of the carriers), the 
formula (|10j) is anisotropic in space in the presence of the spin-orbit coupling, reflecting the symmetries of the crystal 
lattice. If, additionally, the bulk and structure inversion symmetries are broken, it has an antisymmetric part in the 
form of the Dzyaloshinskii-Moriya interactio n 12 ' 13 , 

EE<'E e ---'^'> ( 12 ) 

a" jj' crcr' 

where e^'acr' is the antisymmetric Levi-Civita symbol with e + _ z = 1, while is a pseudovector and exists only in 
the systems with broken inversion symmetry, 



jt — „ fc T"crcr' n-jj' 
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II. SMALL OSCILLATIONS APPROXIMATION 



The system of coupled magnetic moments acts like that of harmonic oscillators, an analogy which is concretised 
mathematically by the Holstein-Primakoff bosonisation^: 

^ « V^l-^Oj , Sr « VSaj^/l-^ , Sj = S - a)a 3 , 

which replaces the spin operators with nonlinear functions of bosonic creation and annihilation operators a}- and a,j . 
To investigate the dynamics of the groundstate and low- lying excitations of the effective Hamiltonian ([9]) , we will use 
the small oscillations approximation and approximate these functions with their power expansions around the state 
of saturation magnetisation: 



5+ w VSa-j , Sj s» VSa}- , Sj — S — a}-a 3 , 



(13) 
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leaving in the Hamiltonian only those terms which are quadratic in the creation and annihilation operators. 

JV 



V 



/km(l /k,m') I z ,2 I t 

"Pkmkm'l | a j a J 



j=l k \ m 



m^m' 



N 



j,j'=l \ kk'mm' 
? 2 JV JV 



/km(l fk'm') i(k'-k)-(R,-R 



V ei (k'-k H R 3 .-R 3 . | 4mk , m ,| 2 (a t ai+a t aj/) 



(14) 



Sfi \^ \^ \^ \^ /km(l /k'm') i(k'-k)-(R 3 — R,Q/- g + 

t/2 2.2.2.^ £ k , m ,_£; km ^ 



' a 'l)^(' S kmk'm' a i + S kmk'm' a '!) 



mk'm' a i + S kmk'm ,u, jV ^kmk'ro' u 3 l, kmk'm'"j"/ 



J — 1 j' — 1 kk' mm' 



Within the small oscillations approximation, the linear and zeroth-order terms do not affect the excitation spectrum. 
This approximation works very well in the long-wave limit, aq -C tt, as the neglected magnon-magnon interactions 
are proportional to (aq) 4 — . 

Firstly, let us take a closer look at the term 



S(3 2 
V 2 



E( /km(l — /k'm') i(k'-k) (R 3 -R,, 

J ,3 —1 \ kk' mm' 



^l s kmk'm'| 2 ) a j a j 



Summation over j' gives X^=i e ^ k k - ) R ' i ' = -/V<5kk', which allows us to write the above term 



as 



NS/3 2 
~V 2 ~ 



N 



£ EE 



/km (1 /k-rn' ) | 
' -Ek,m' — ^k 



5 kmkm'| 2 I a j a 3 



j—l \ k mm' 

The twin term with aLa/ at the end is treated in the same fashion, after which we write ([140 as 



j — 1 k \ m m^m 



N 



/km(l /k,m') | 
, ^k,m' — ^k 



s kmkm'l I a j a j 



2NSf3 2 



N 



E EE 

j — 1 \ k mm' 

, 2 JV JV 



r 2 



/km(l /k,m') i z 
— p ^""p l S kmkm'l I u j u J 



g/3 " " ^ ^ / km (l - /k'm') j(k'-k)-(R,-R,.,)/ + «. + .- at)tf s + a . +s " a t) 

t/2 2-^2-^2-^2-^1 TP —TP \ b Wmk'm' u 3 < °kmk'm' a j) l\mk'm' "j ^ s kmk'm' U J > ■ 

v -1,111/ , -^k'm.' -t^k™ 

j — 1 j' — l kk mm' 



Due to the equality A = NSfi/V, terms with |s k mkm'l 2 cance l an d we get 



JY 



^ ff = -7(EE«^]E^ 



k m 



y2 2^ 2^ 



/km(l fk'm') 



kk' mm' 



N 



E el< ' k k ^ Rj ( S kmk'm' ffl J + S kmk'm' ft j) 



3'=1 



JV 



E £l< ' k k ^ R/ ( S kmk'm' a i' + s kmk'm' a ]') 



We can now proceed with the Fourier transform (invoking the virtual-crystal approximation), using the fact that 

EN f t 



^° ff — ~T7 \ 2\2i E / kmS kmki 



V 



k m 



, E fl q' 

/ q 



/km(l fk'm') 



S kmk'm' ak '- k + S kmk' 



kk' 



J k'-k) ( S kr) 



,k'm' ak '- k + S kmk 



'm' u k'-k y 
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After simple algebraic transformations, including the symmetrisation of the sums over wavevectors, we arrive at 
the final form of the harmonic Hamiltonian, 

"W cff = Xl[( S ™ *q~) a q fl q- ^Xq + a q a- q - ^X q _ a q a -q • ( 15 ) 
q 

We call it the interaction representation as it describes the perturbation of the ground state by the isotropic Coulomb 
interaction (first term) and by the spin-orbit interaction, coupling modes of different q (remaining terms). The spin 
susceptibility of the carriers is given by 

,,<x<x' _ nS P 2 ST ST f km ~ / k +q. m ' a <r c cr' /1a s 

Aq — T/ 2^1 VP _ P A km(k+q)m' ; '(k+q)m'km ) (iDj 

k mm' ^ 

where n = N/V is the density of localised spins S in the sample volume V and nS/3 = A. The presence of the energy 
denominator shows that Xq CT corresponds to the second-order part of the Hamiltonian ([7]). As promised in Sec.|H the 
vanishing of the denominator is not harmful, due to the de l'Hospital rule. 

The formula (fT5|) implies that x+ + = (x q ~)* is symmetric in q. In the absence of the spin-orbit coupling the 
bands' spins become fully polarised, which causes Xq + to vanish, as (4>\s + (ip' \s + \iJj) = for any choice of spinors 
ip,ip' = |t), This is also true for non-zero spin-orbit coupling in the case when the valence bands are isotropic 16 
and for s-type bands in general, as for them the total angular momentum is equal to spin and thus they are fully 
polarised even in the presence of spin-orbit interaction. 

Because x+ ~, x q + S ^ inherit the symmetry of the ipkm eigenstates, it can be expected to be symmetric with 
respect to q for systems which preserve space inversion symmetry, like in the case analysed in Ref. [l6l and otherwise 
for systems which do not. The q- independent term describes the interaction of a single magnetic ion with a molecular 
field arising from the intraband spin polarisation of the carriers, 

a = — y /k,m 4mkm ■ (1?) 

k m 

The corresponding term reflecting the interband polarisation, 

_ nSp \^ \ ^ /km ~ /k,m' | a ,2 /,o\ 

"so— Tj— 7j 7_j 1? 5 l*kmkm'l > \ L0 ) 

V , '-T . &k.m' — fikm 
k mfm 

arises from both Hjp (|TU)) and the second part of Hj coefficient (TTT|) . and cancels exactly in the full Hamiltonian 

"H cff . As announced in the previous section, this shows that the Bloembergen-Rowland exchange, driven by virtual 
spin transitions between different bands, does not exist in the low-temperature limit covered by small oscillations 
approximation. 



III. BOGOLIUBOV TRANSFORM 



So far, we have used the Lowdin calculus to find the effective Hamiltonian of lattice ions interacting through 
delocalised carriers in a self-consistent way, written in the interaction representation of creation and annihilation 
operators. To calculate the dispersion dependence of its low lying energy states, spin waves, we need to transform the 
Hamiltonian to the spin-wave representation. For this purpose we will use the Bogoliubov transform, which we have 
adapted to systems with inversion symmetry breaking. 

The effective Hamiltonian for the lattice ions derived in Sees. UlandlHlin the interaction picture has the following 
form: 

n ° S = H [( S ~ *q~) a q°q " 2 X ^ +aqfl ~ q ~ 2 Xq ^ aq<I - q ' 
q 

It describes the spin system in terms of circularly polarised plane waves (the first term), which interact with each 
other and deform in time (the remaining terms). 

We want to obtain the dispersion relation of independent, stable magnons. For this purpose, we diagonalise H cS by 
the Bogoliubov transformation from a q , a q to 6 q , & q operators (which describe independent excitation modes), keeping 
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in mind that we deal with the system which breaks the space inversion symmetry. Because 'H off in the above form 
mixes states with opposite wavevectors, we write the sum over q in an explicitly symmetrised form: 



H ° S 2 

q 



\ E [( S ~ X q ) a q°q + ( S ~ X-q) a -q a -q - 7j(Xq + + xt q )a q a_ q - ^(x q + X- q )a q a -q 



(19) 



The Bogoliubov transformation is given by the formula 

fl q = U q fe q + W q 6L q , U q , »q £ C , (20) 

with the conditions 

1 = |Uq| 2 - |Uq| 2 , = U q V- q - U q 1i_ q (21) 

ensuring the preservation of canonical commutation relations (equivalently, the invertibility of this transformation). 
Inserting this into (|19[) . we obtain (neglecting the lower-order terms and making use of the properties of x's) 

q 

+ ( S - X-q) (|«-q| 2 6tq&-q + l^-ql^q^q + "-qf-q^qftlq + U-qU* q 6 q 6_ q ) ( 2 2) 

- Xq + (w q U_ q 6 q 6_ q + Uq«-q&q&q + fqW-q^ -q&-q + "qW-q&q&lq) 

~ (Xq + )*( U q U lq^q^-q + «q«tq & q & q + W q U lq^ q &-q + W q U*q&q&-q) , 

which leads to the following complex equation for each q: 

(S " Xq~) «q«q + (S " X-q) «-q«-q " Xq + «q«-q " (X + + )X«-q = . (23) 

Since all equations for zt's and u's are invariant under the reflection of q, we assume that u q = u_ q and w q = v_ q . 
The standard parametrisation for u q , w q consistent with the first condition (|21[) reads 

M q = e lMq cosh q , w q = e Wq sinh q . 

Inserting it into (|2"5j) gives 

0= (2S-Xq _ -Xt q ) cosh q sinh 0qe i O**-"'O C osh 2 q e 2 ^ - (x+ + )* sinh 2 6»qe~ 2 *^ 

= (2S - x q ~ - Xtq) cosh q sinh q - x q + cosh 2 q e i( ^ +,y - ) - (x+ + )* sinh 2 q e- l ^ + ^ . 

There is no equation for /z q — so we assume it is zero. We can replace /x q by another unknown, /i q , defined by 

Xq + e 2i ^ = IXqV^ , 

resulting in 

(2H - xt ~ X^) cosh0 q sinh0 q = |x+ + | (cosh 2 q e 2 < + sinh 2 q e- 2 <) 

or 

(2S — x q ~ — X-q) cosh q sinh q = |x q + l (cosh 2 q + sinh 2 q ) cos2/i q + isin2^ q , 
which requires fx' = 0, leading to the equation 

(2S - Xq~ - X-q) COSh0 q Sillh q = |x+ + |(cOSh 2 q + Sillh 2 q ) , 

hence 

20 q = arctanh — , (24) 
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with the solution non-existent when 



lx+ + l 



> 1. This particular- - solution will turn out to be sufficient for 



our needs, as numerical calculations give |Xq + | ^ I s — s(Xq + X-q)|> 
Gathering all terms multiplying 6 q 6 q in (|22|) . we obtain 



1 r 
2 
l 

2 



q i 2 + (s-x+ q ; 



Kql 2 - Xq + M q«-q ~ 0Cq + 



5 - x+-) cosh 2 q + (S - X + q ) sinh 2 q - 2| X ++| cosh# q sinh0 q 



„ Xq +X- 



-| Xq + |sinh20 q 



Using (|24p and assuming that H - q 2 q > (which is the case for our numerical calculations), we obtain the final 
form of the effective Hamiltonian in the spin-wave representation: 



q > 



(25) 



where excitation modes are spin waves with dispersion 19 

X-q — Xq 



/(2S- Xq — X- + q ) 2 



(26) 



In the case of Xq = X-qj fulfilled for the systems invariant under space inversion, the above formula simplifies to the 
solution by Konig et al~, 



(S-Xq- 



IXqT 



Furthermore, neglecting the spin-orbit coupling, when x ++ = 0, Bogoliubov transformation is unnecessary, and the 
effective Hamiltonian is already diagonalised by a q , a q operators. 
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In general, P cx is not a projection operator, as (P p + |P S ) 2 = P 2 + f(P p P s + P S P P ) + frP? =Pp+ jpPs + P cx for a / ft. 
However, it maps its image onto itself. 

Other solutions may exist which do not fulfil our assumption that u q = u_ q and i> q = ii_ q . However, for a finite number 

of degrees of freedom the representation of canonical commutation relations is unique up to a unitary transformation, 

which guarantees the uniqueness of the obtained dispersion relation. In general, a set of creation and annihilation operators 

diagonalising H eS must exist, but it does not need to be obtainable via a transformation of the form p0[) . 

The negative sign of the Hamiltonian (|25l) is a consequence of using the electronic convention to describe the carrier-ion 

system. 



